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Abstract 

A method of solving Maxwell equations in a vicinity of a multi- 
pole particle (moving along an arbitrary trajectory) is proposed. The 
method is based on a geometric construction of a trajectory-adapted 
coordinate system, which simplifies considerably the equations. The 
solution is given in terms of a series, where a new family of special 
functions arises in a natural way. Singular behaviour of the field near 
to the particle may be analyzed this way up to an arbitrary order. Ap- 
plication to the self-interaction problems in classical electrodynamics 
is discussed. 



1 Introduction 

Classical approach to the problem of motion in classical electrodynamics is 
due to Dirac (cf. [3], [6]). In this approach, equations of motion of charged 
point particles, interacting via Maxwell field, are derived from field equations 
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by the following procedure: a generic state of the composed "particles + field" 
system is treated as a deformation of the "ground state" of the field, uniquely 
determined by the positions of the particles. The ground state is defined 
via a (non-local in time) decomposition of the actual field into the retarded 
(or advanced) field and the remaining "radiation field". Unfortunately, to 
decide what is the retarded or the advanced field, the entire trajectory of the 
particle must be known in advance and, whence, the causality of the theory 
is violated. 

Recently, it was shown that such a non-causal procedure may be avoided. 
In this approach, equations of motion are a simple consequence of the conser- 
vation laws imposed on an (appropriately defined, "already renormalized" ) 
total four- momentum of the "particles + field" system (see e.g. [7], [5], [TT]). 
Here, the "ground state" of the system may be defined via a conditional min- 
imization of the energy, with the positions and the charges of the particles 
being fixed. 

Mathematically, this leads to a simple (elliptic) variational problem for 
the behaviour of the field in a topologically non-trivial region of K 3 (exterior 
of the particles), where the charges of the particles provide the necessary 
boundary conditions. Such a reduction of electrodynamics to electrostatics 
provides basic information about the behaviour of the field in the vicinity 
of the particle (e. g. the leading r~ 2 term for the monopole particle or the 
r~ 3 behaviour for the dipole) but fails to capture its more subtle features, 
which are necessary to describe such phenomena like polarization, which is 
necessary for the analysis of the dynamical stability of the classical elec- 
trodynamics (see e. g. discussion of these issues in [10J and the references 
therein). 

A substantial improvement of description of the field singularity near 
the trajectory is obtained if we reduce the dynamics with respect to a one- 
parameter group of boosts (see pj] ) , instead of the one-parameter group of 
time translations, corresponding to the standard, flat "(3+l)-decomposition" 
of the Minkowski spacetime. As a result, we obtain the "electrostatics in 
a uniformly accelerated reference frame" as a proper tool to construct the 
ground state of the field. Here, field equations are again elliptic, but the 
flat Eucledean geometry is replaced by the Lobaczewski geometry. Within 
this framework, the Born solution describing Maxwell field attached to a 
uniformly accelerated particle (see e. g. [13]) arises as a fundamental solution 
of the Laplace-Beltrami operator in the Lobaczewski space: A Lip = So, 
where the right hand side is the Dirac delta distribution (a nice analysis of 
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this structure was proposed by Turakulov in [H]). 

In the present paper we show how to use this (naturally arising) Loba- 
czewski three-geometry to the analysis of the full dynamical problem: analy- 
sis of the Maxwell field generated by a (moving along an arbitrary trajectory) 
multipole particle. Our method consists in splitting the portion of spacetime 
in the vicinity of a generic trajectory of a particle into a collection of 3- 
spaces orthogonal to the trajectory Each of these spaces carries a natural 
Lobaczewski geometry, uniquely implied by the instantaneous acceleration 
of the particle. As a result of this construction, the flat spacetime geome- 
try is replaced by a curved (general relativistic) geometry, whose parameters 
are uniquely defined by the parameters of the trajectory in question. We 
show that the Maxwell equations may be rewritten in terms of this geome- 
try. This method provides a new, powerful tool, which enables us to solve 
Maxwell equations by a successive approximation method. 

The results presented here give us an important improvement and gen- 
eralization of earlier methods (see [8] or [9]), based on the use of the Fermi 
coordinates, defined by the Fermi tetrad attached to the trajectory. The 
geometric construction proposed in the present paper will be referred to as 
a modified Fermi system. 

Naively speaking, the standard (3+l)-approach consists in approximating 
the generic trajectory by a straight line, matching only the instantaneous 
velocity of the particle. In the approach based on Fermi coordinates, we 
use much better approximation given by hyperbolic (uniformly accelerated) 
trajectories, which match not only the velocity but also the acceleration of 
the particle. The value of acceleration is encoded in the parameter of the 
Lobaczewski 3-space arising in this construction. In our new approach (which 
we refer to as the "modified Fermi"), proposed in the present paper, we use 
still better approximation: instantaneous acceleration a(t) of the trajectory 
implies the parameter of each of the (Lobaczewski) 3-space folia {t =const. } 
of the constructed geometry; there is, however, a non vanishing shift vector 
whose value encodes the derivatives a^(t). 

The basic technical tool of our method consists in deriving an explicit 
formula for the operator inverse to the Lobaczewski curl operator (which 
we denote by "curLj") on a certain space of divergence-free vector fields. 
The inversion requires a successive solution of Legendre differential equation 
with non trivial right hand side. This, in principle, could lead very quickly 
to non-elementary functions (see formula fl33|) in Section H]). To our great 
surprise, we were able to express the results (up to the order in the 
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radius r) in terms of rational functions of two universal combinations of the 
radial coordinate, namely u = ar/2 and R = log(^pjj), where a is the (scalar) 
acceleration of the particle (see Section [6]). This result (together with the 
extremely simple and natural physical context) leads us to a conjecture that 
the functions obtained this way form a new family of special functions, which 
we discuss at the end of the paper. 

Physically, our results form a basis for the theory of motion of a po- 
larizable particle: a particle which carries a dynamical electric or magnetic 
moment, i. e. a moment which is not frozen within the particle but depends 
dynamically upon the surrounding electromagnetic field. There is a serious 
hope that such a theory would be free of the standard non-stability prob- 
lems accompanying the Abraham-Dirac theory (see [10J and the discussion 
therein). 

The paper is organized as follows. In Section [2] we define the modified 
Fermi system. In Section 3 we derive Maxwell equations in the modified 
Fermi system and interpret them in terms of the Lobaczewski geometry. In 
Section 4 we define basic fields X and Y and show that the operations div, 
curLp and # act on X, Y in a simple way. These results give a generalization 
of the construction proposed in [9] . In Section 5 we solve Maxwell equations 
in a vicinity of the multipole particle, in terms of a series consisting of the X 
and Y fields. The procedure is analogous to that used in [9]. In case of the 
monopole and the dipole particle, the resulting series up to the order r 1 is 
explicitely given in Section 6, where r is the radial coordinate in the modified 
Fermi system. In Section 7 we propose a conjecture concerning the special 
functions arising naturally in our construction. Several specific computations 
and proofs are shifted to the Appendices A-D. The results of Sections 6 and 
7 were obtained partially with help of the symbolic computations program 
MAPLE 8. 

Thorough the paper we use greek indices (running from to 3) to label 
space-time coordinates, and latin indices (runing from 1 to 3) to label space 
coordinates. We always use Einstein convention: summation over repeated 
indices (in both space-time and space cases). The components of an n-tensor 
T are denoted by T ili2 „ An . 
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2 The Fermi— propagated and the modified 
Fermi— propagated systems 

In this section we recall properties of the Fermi-propagated system of coor- 
dinates, in which the particle "remains at rest" at each instant of time. Next 
we modify that system in a way which is similar to transition from cartesian 
to bispherical (bipolar) coordinates in M 3 , with the particle's position being 
one of its centers (poles). 

Let y x , A = 0,1,2,3, denote the (Minkowski) spacetime coordinates in 
a fixed inertial ("laboratory") system, corresponding to the metric tensor 
77 = diag(— , +, +, +). Consider an arbitrary particle's trajectory q x (t) = 
(t,q k (t)). Let t = r(t) denote a particle's proper time along the trajectory. 
The normalized four- velocity is given by u = dqjdr and the particle's accel- 
eration by a = du/dr = d 2 q/dr 2 . We define (see e.g. [11]) the rest-frame 
space £ T ( t ) as the hyperplane orthogonal to the trajectory (i.e. to the vector 
e(o) := u) at the point q(t). We choose an orthonormal basis e^, I = 1, 2, 3, 
in S T , such that e^) are positively oriented. Thus (e^\e^) = rj a p. Since 
(u|a) = 0, one gets a = a l e^ for some a 1 . 

At each instant of time r, the above system may be subject to the SO (3) 
group of rotations of the Dreibein (e^^r)), playing role of the group of gauge 
transformations. The Fermi-propagated system (see e.g. [UJ) is defined by 
the following condition: 

de^/dr = a l u . 

Once the above condition is satisfied, there remains a single, global (time- 
independent) SO(3) gauge transformation, corresponding to an arbitrary 
choice of initial conditions: e^To). Moreover, we obtain in this case: 

du/dr = a ern , 

which enables us to interpret a 1 , I = 1,2,3; (with a = 0) as components of 
the acceleration vector. Next, we define in a neighbourhood of the trajectory 
the Fermi-propagated (local) coordinates £ M := (£°,£ z ) putting £° := r and 

y X = q X (r)+ee^(r). (1) 

Hence, are cartesian coordinates on S T , related to em and centered at the 
particle's position (£ = for y = q(t)). The Minkowski metric tensor 7 in 
that system is given by 

Iki = hi, loi = 0, 700 = -N 2 , (2) 
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where k, I = 1, 2, 3, N = 1 + a/£' (see e.g. [TT], p. 373, but with £ denoted 
there by x). 

Now, we change variables on each S r and define the "modified Fermi- 
propagated" system x M = x l ) by putting x° := r = £° and: 

x l = ^ + ¥(? (3) 

where a = (jki^a 1 ) 1 ^ 2 , p = ('jki^F, 1 ) 1 ^ 2 - It is easy to check that the inverse 
transformation is given by the following formula: 

/ 1/2 

j x — !<rr z 

£ 1 i I 1 9 9 ' 

1 — a,j£ l + ^aT^ 

where r = (7/c/x fc x / ) 1 / 2 (a, = a* is the same as before because the particle is 
situated at x = £ = 0, where dx k /d£} = 5 k and (j2J) holds). 

A simple computation shows that the flat (Minkowski) metric tensor has 
now the following components: 

N 2 d£ k d£ k 2 di k di k 

9kl — —fOkl, gol — ~5 "fTT; #00 — — ^ + ~^ q > (5) 



(p z or ox 1 or or 



with 



(p = 1 - (ar/2) 2 (6) 

Of course, we have: g 00 = g(dr ; dr) = —N~ 2 , and thus the lapse function N 
is the same as for the previous Fermi system, because the time variable has 
not been changed (cf. [9]). In terms of the new coordinates, its value may 



be written as follows: N = (p/M, where we denote M = 1 — aix 1 + ^a 2 r 2 
(1 + a k £ k + \a 2 p 2 Y x . Moreover, y/det g kl = (N/pf = M~ 3 . Denoting by 
g the 3-dimensional inverse of (gki), k,l = 1,2,3, the shift vector may be 
calculated as follows: (cf. [9]) 

N k = ? l g* = |^ = - \r 2 d k , (7) 

where "dot" denotes the derivative w.r.t. r, e.g. dj = daj/dr. 

To illustrate the geometric structure of the above coordinate system, let 
us apply to the above picture such a rotation O, which positions the third axis 
e(3) in the direction of the acceleration a, i. e. such that we have: a = a e^)(r). 
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Apply now rotation O to coordinates (x 1 ,^ 2 ,^ 3 ) and denote by (z 1 , z 2 , z 3 ) 
new coordinates on E T , obtained via such a rotation. Next, using z k as Carte- 
sian coordinates, construct the corresponding spherical coordinates (r,r],(f)). 
Finally, define the variable \i by 

r = - exp(— u). 
a 

Then, it is easy to show that (/i, r], <p) form the bispherical system of coor- 
dinates on the Euclidean space £ r , with the particle's position being one 
of its centers: r = (i.e. /i — > +00). The vector connecting this center 
with the other center of the bispherical system is parallel but opposite to the 
acceleration and its length equals -. 



3 Maxwell equations 

Maxwell equations: df = and d*f = J, with * denoting the Hodge 
"star" operator, can be written in an arbitrary system of coordinates, in an 
arbitrarily curved spacetime, as follows (cf. e.g. [9J): 

[7 / H = 0, (8) 
d v T^ = J», (9) 

J*» = V^teTg sT^faft , (10) 

where the brackets "[ ]" denote the complete anti-symmetrization. More- 
over, we denote: = y/— det g J M , where J 7 ' is the four-current vector 
and J 11 is the current density (an "odd three- form", see e. g. [2]), sat- 
isfying the continuity equation d^J^ = 0. The two vector-densities: V 
(electric induction) and B (magnetic induction), are defined (cf. [9]) as the 
following components of these tensors: \fki^ klm — B m , jF ofc = V k , where 
e kim j g g^andard Levi-Civitta symbol (totally antisymmetric and nor- 
malized: e 123 = 1). Due to (jSJ), the covariant components of the mag- 
netic field, calculated in our modified- Fermi coordinates (ED, are equal to 
(cf. [9j): B m = g mk B k = g mk (det g kl )- 1 ' 2 B k = { V /N)B m . Hence, we have: 
NB m = ipB m and, similarly, ND m = (pT> m . These formulae show that a sub- 
stantial simplification of the field equations is obtained, if we use a fictitious, 
flat metric 8 k i to identify "upper and lower indices", i. e. to fix an isomor- 
phism between vectors and covectors. From now on, we strictly observe this 
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convention. Consequently, we can rewrite Maxwell equations (151)- (ITU1) (cf. 
(9)-(12) of 0) as equations for the two quantities (V,B), treated as two 



vectors in this fictitious flat geometry: 

d k V k = J\ (11) 

d k B k = 0, (12) 

V k - di{N l V k - N k V l ) = e m d i {<p8i j B j ) - J k , (13) 

B k - d t (N l B k - N k B l ) = -e kU di((p5ijV3), (14) 



where the vector N k is uniquely implied by the derivative of the trajectory up 
to the third order, according to formula (171), whereas the "conformal factor" ip 
is implied by the second derivative, according to formula ([6]). We introduce 
the following short-hand notation for the differential operators appearing 
here: 

W # =W + W, (15) 

where 

W k = di(N k W l - N l W k ) , (16) 

and 

curl^W := curl(^W) . 
This enables us to rewrite Maxwell equations as follows: 

divX? = J°, (17) 

divB = 0, (18) 

V* = curing - J, (19) 

E* = -curl^V . (20) 

Now, we are going to show that the differential operator "curl^" can be 
nicely interpreted as a genuine curl operator in the Lobaczewski geometry. 
For this purpose consider the covectors which are subject to differentiation 
in the definition of curl^,, see formulae (fT3"]) . f[T4l : 

d { := (fSijV 3 , hi := <p5ijB j . 

We want to find an appropriate (yet another) 3-metric \x on each leaf E r , 
such that the above covectors (differential one-forms) are related with the 
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vector densities D and B (differential three-forms) via the standard Hodge 
formula: 

V k = y/Jl^dt , B k = y/fi^bt . (21) 

It is easily seen that /i must be conformally similar to the flat metric: = 
XSki- Taken into account that ■J]ifJ l kl = X^8 kl , we conclude that to fulfil (f2~TI) 
we must take cpX? = 1, i. e. A = (p~ 2 . It is a matter of easy calculations to 
check that, indeed, the resulting metric tensor: 

Hid : = (f~ 2 S k i = - / 5 H , (22) 

1 — {ar/2) z 

defined on the ball K := {x G IR 3 : r = < -} is equal to the Lobaczewski 
metric with constant negative curvature R = —6a 2 . 

In the particular case of a uniformly accelerated trajectory, i. e. when 
a — 0, the shift vector vanishes: N = 0. As noticed by several authors, 
the "electrostatics" (with respect to such a uniformly accelerated reference 
system!) implies: curl^P = 0, equivalent to d — grad ip. Then, Gauss 
equation (|17p reduces such a theory to the standard "potential theory" in the 
Lobaczewski space: div grad = A^ip = J° . The fundamental solution of 
this equation, corresponding to the point particle: J Q = 5 (Dirac-delta-like 
charge density - cf. [H]), describes the electromagnetic field accompanying 
the uniformly accelerated particle and is equal to the Born solution (cf. [TJ 
or P3]) 



4 The calculus of X and Y fields 

In this section we define elementary fields X and Y, which enable us to 
reduce the asymptotic analysis of Maxwell equations in a neighbourhood of 
a moving, polarized particle, to a relatively simple algebra of these fields. 

Given a pair (n, m) of natural numbers such that n = 1,2,..., m = 
0, 1,2,.. ., consider the vector space L n m of Laurent series of a single variable 
r, of the form 

oo 

Y^c k r' n+2m+2k , 

k=0 

convergent in an annular neighbourhood < r < e, where the coefficients c k 
are real. Observe that we have L n m+1 C L n ,m- By the order of the series we 
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understand the order of its first non-vanishing term (i. e. / is of order I if it 
behaves at r = like r l ). 

For any pair (/, Q), where / G L n0 , and Q = (Qi r ..i n ) is a completely 
symmetric, traceless tensor of rank n (i. e. satisfying: 6 %in Qi 1 i 2 ...i n = 0), 
consider the following vector fields: 

Xf n) (/, Q) = [(n + 1)/ - r/, r ] ^x k Q h ... in x^ •••**" + f )T L Q l. nX n . . . a A. j 

(23) 

Y { i ) (fM) = ^^ ri e kll x,Q U2 ... tn x^---x^ . (24) 

For an exceptional case n — 0, every Q is a scalar and, therefore, the field 
Y(o) is not defined. Moreover, we assume in this case that / is constant. 

This is an appropriate generalization of the fields X and Y, which were 
defined in [5] for purposes of the analysis of the Maxwell equations in terms of 
the standard Fermi system. As will be seen in the sequel, this generalization 
enables us to describe the asymptotic behaviour of the Maxwell field in a 
neighbourhood of a polarized, point particle in a much efficient way. 

It can be easily shown (cf. Appendix A) that, outside of the center r = 0, 
the fields X and Y are divergence-free: 

divX (n) (/, Q) = divY (n) (/, Q) = 0, (25) 

and the operator curl^, acts in a simple way on these fields: 

curl„Y (n) (/, Q) = -X (n) (/, Q) , (26) 

curl^X (n) (/,Q) = -a 2 Y (n) (V)(/),Q), (27) 

where 

h(n){f) ~ —tft,rr + y/,r + ^ / ■ (28) 

(The full action of these operators understood in the sense of distributions 
contains also Dirac-5-like terms at r = 0, cf. Appendix [D]) . A much simpler 
form of the operator h( n ) is obtained if we introduce the following new vari- 
ables: u = ar/2, z = \{u~ l + u) — + ^) = cosh/i. Then z > 1 and as- 
suming u < 1 (i. e. near to the particle) we have: r = 2u/a = -(z — \/ z 2 — 1). 
Now, h( n ) takes the form: 

h (n) [g(z)} = (z 2 - 1){(1 - z 2 )^ - 2z^ + n(n + 1)^}. (29) 
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Equation h( n )[g(z)} = is equivalent to the Legendre equation and, therefore, 
has two independent solutions: the Legendre polynomial of order n: P n (z) = 
lh. S^ 2-1 )™ = A ^ r ) and the Legendre function \P n (z) log(f^)+w n _i(z) = 
B n (r), where v n -\ is a polynomial of order n — 1, V-\ = 0, P n and v n are even 
(odd) for n even (odd) (cf. [13], pp. 302 and 319). Thus A n (r) and B n {r) 
are solutions of equation h^(f) = 0. A n has the same parity as n, B n - the 
opposite parity (one can formally define them for \u\ — 4p < 1). 

It can be proved (see Appendix [B]) that A n is a Laurent polynomial of 
r with the lowest power r _n , whereas B n is regular at r = and its Taylor 
series starts from r n+1 . 

Proposition. There exists a unique operator (the index m is suppressed 
for simplicity) s (n) : L nm -> L n)OT+1 , which is inverse to h {n) : L„ im+1 -> L nm . 

Proof. Because solutions of the homogeneous problem are known, equa- 
tion 

h(n)(f) = l (30) 
can be solved by the "variation of constants" method: 

f(r) -A(r)f a %)^(r)dr _ f a 2 l(r)A n (r)dr 

fir) - A n{ r) j ^ 2(Aai _ B n{ r) J _ (31) 

where the Wronskian A n B n ' — B n A n ' is proportional to y? -1 , so the denom- 
inators under integrals are proportional to cp. Assume now that / 6 L nm . 
We define sr n )(l) as the right hand side of fl3~T|) with the integrals performed 
term by term, using the formula J r a dr = (a ^ —1). Then S( n ) is linear. 
Moreover, the Laurent series of the function under the first integral is odd 
and starts at least from r 1+2m , so the series for the integral is even and starts 
at least from r 2+2m . The Laurent series under the second integral is even and 
starts at least from r ~ 2n + 2m . its integral is odd and starts from r - 2n + 2m + 1 . 
(This proves that, indeed, the case a = — 1 never occurs.) Thus, S( n )(0 is 
even or odd (depending upon n) and starts at least from r _ n+2 ( m +i\ i. e. we 
have S(n){l) G L n>m+ i. Any other solution / of ( 1301 differs from S( n )(0 by a 
combination of A n and B n . But the condition / G L n ,m+i excludes both A n 
(its series starts from lower power r _n ) and B n (its parity is different from 
n). This proves the uniqueness of S( n ). Q.E.D. 

We conclude that, restricted to the space of those X( n )(/, Q) and Y( n )(/, Q), 
for which / G L nm , the operator "curl^" has a right inverse "curl" 1 " (i.e. 
curl^curl" 1 = id) given by the following formulae: 

curl; 1 X (n) (/,Q) = -Y (n) (/,Q), (32) 



11 



curl^Y^/, Q) = -a- 2 X (n) ( S(n) (/), Q) . (33) 

Let us notice that for / e L nm being of minimal order (— n + 2m), one has 
that: 

X (n )(/, Q) is at least of order r~ n+2m - 2 , 



Y (n) (/, Q) is at least of order r' n+2m -\ 



\n){ 

Those orders are called the generic orders. Observe that (generically) curl" 1 
increases the order (in r) by 1. 

Now, suppose that both / and Q are time-dependent. It turns out that 
also time-depending operator # (cf. (TT51) - (IT6"1) ) acts in a simple way on the 
fields X and Y. Indeed, we have: 

Theorem 1. For r ^ 

X W (/, Q) # = -^X (n+1) (r 2 /,r + (n + l)rf, b V Q) + Y {n) (g h b x Q) + 
n 2 - 1 

+ 2n(2n + 1) X (n _ 1) (-r 2 /,r + nrf, bJQ) + X (ft) (/, Q) + X (n) (/, Q),(34) 

Y (n) (/,Q) # = Y (n+1) ( 92l bVQ) - - — X w (^,bx Q) + 

2{n + 1) <p 

+Y (n _ 1) (< 7 3, bJQ) + Y (n) (/, Q) + Y (n) (/, Q) + Y (n) (^, adQ), (35) 

where bk = (ofe)", a = da/dr = a h bk/a, whereas the symmetric traceless 
tensors: b V Q, b x Q and bJQ of ranks (respectively) n+ 1, n and n — 1, 
are built from the n-tensor Q and the vector Q according to the following 
formulae (the letters placed below dots denote missing indices): 



1 n 

(b V Q)i ...i n = — — V b is Q io in 
n + 1 

s 

2 m 
(n + l)(2n + 1) ^ W> Q«*o^i», (36) 

jw 

1 n 

(b x Q) h ... in = - £ e^ m & ; Q mil (37) 

k 

12 



9i 



(bjQ)j 2 ...j n b Qmi2---i n i 
<P r 2. 



(38) 



2(n+l) 



[r7.rr + 2r/, r -n(n+l)/], 



92 = — (ft + 3)r f r f r . 



n 



93 



n(4n + 2) 



(ft - 2)r/ - r 2 /, r - 



a 2 fr 3 
2<p 



(p = l- (ar/2) 2 . 

Remark. For the properties ofbVQ, bxQ, bJQ (computed in a simpler 
situation, but valid universally) see [9], pp. 304-305 (Warning: Normalization 
of the tensors b V Q and b x Q used in fl36l) - (l37|) differs slightly from that 
used in (23)-(24) of [9J). Let us notice that for W being X or Y field, W # 
is at least of the same generic order of r as W. If, moreover, /, Q do not 
depend upon r, the generic order of W* increases by at least 1 with respect 
to the generic order of W. 

We shall use these formulae in a special case, when / depends upon time 
only via the acceleration a = a(r) contained in the combination u = ar/2. 
Assuming, therefore, that f(r, r) = }'{u), one can rewrite (I3"4"l) - (13"51) as follows: 



X( n )(/, Q) # = -X( n+ l) 



u 2 f, u + (n + l)uf, -b V Q 



+ 



———Y {n) [ V {u 2 f\ uu + 2uf\ u - n{n + l)/),b x Q] + 



n 2 -l 
n(2n + 1 



X 



(n-l) 



-u 2 f >u + nuf, -bj Q 
a 



+ X (n) [ U /, tt ,-Q]+X (n) [/,Q],(39) 
a 



Y (n) (/,Q) # = -Y 



(n+l) 



(n + 3)uf + u 2 f, u + 



2u 3 f 1 



-bVQ 



+ 



n + l 



X 



(») 



n 2 - 1 
n(2n + 1 



Y 



(n-l) 



(n - 2)uf - M 2 /,n 



(p a 
,fu 2 1 

-,- bx Q + 

<p a z 
2u 3 f 1 



-bJQ + 

(p a 

a • fv 2 2a 

Y(n)[uf, u , -Q] + Y (n) [/, Q] + Y (n) [^— , — Q]. 

a ip a 



(40) 



The proof of Theorem 1 is given in Appendix A. 
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5 Expansion of the Maxwell field 



In this section we find an asymptotic solution for the Maxwell field in a vicin- 
ity of a time-dependent, electric 2 n — pole, moving along a given trajectory. 

Given a timelike trajectory, we construct its modified Fermi frame like in 
Section 2 and consider Maxwell equations (TITll - lpOl) with electric charge: 

J° = c n Q i ^d il ...d i J^ (41) 

and the corresponding electric current 

J k = -c n Q ki *- i "d l2 ...d l jW, (42) 

at the right-hand side. The above four-current is conserved: d^J^ = 0, cf. 
[9]. Here, c n = (— l) n AiT/(2n — 1)!! is the normalization factor. For n = 
we assume J k = and Q = const. (See Appendix C for relation of these 
2 n -poles, defined in the modified Fermi frame, with the standard 2 n -poles in 
the Fermi frame.) 



Theorem 2. Maxwell equations flUD-tEDD witn the sources (jH])-(jI2]) are 
solved by the following (formal) series: 

V = V d + P ( _ n _ 2) + £>(_ n) + X>(-«+2) + • • • , (43) 

B = + B ( _ n+1) + B ( _ n+3) + . . . , (44) 

where 

K = (-l) n (2n + 1)! , ^ 2 -"^ • • • d*J (3 \ (45) 

P ( _„_ 2) =X (n) (k n P n {z),Q), (46) 
% + i) =curl; 1 [P (fc) # ], k = -n-2,-n,-n + 2,..., (47) 
%•+!) = -curl- 1 [i3 (fc) # ], fc = -n-l,-n + l,-n + 3,..., (48) 

i.e. B(_ n _i) = curl- 1 [D ( _ n _ 2 ) # ], £>(-») = -curl" 1 [£(_„_!)#], 

S(_rH-i) = curl^ 1 [X> ( _ n ) # ], fc n = (2a) n /( 2 ™). Moreover, B ( q and £> (i) are 

at least of order r l . 

Remark 1. The fields B, T> and the operations div, curl^, # are consid- 
ered here in the sense of distributions. On the other hand, curl" 1 acts on 
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the fields X and Y via (I521)- (I55|1 . One can add to flT7j) - (120]) any regular 
solution (V, B) of homogeneous Maxwell equations. 

Remark 2. One has k n P n (z) = r~ n + /„, where lim r ^o f n (r)r n = 0. 
Thus the leading term in V equals T>d + X( n )(r~ n , Q) = — grad^n where 



is the standard potential for the classical 2 n -pole ( l4"Tj) with a = 0, Q = const, 
(cf. Appendix D). 

Proof of Theorem 2 is given in Appendix D. 

6 Application to the monopole and the dipole 

cases 

In this section we apply the above formulae to the case of a monopole and a 
dipole particle. We find explicitly asymptotic expansion of the Maxwell field 
(T>, B) up to the order r 1 . 

In the monopole case we start from £>(_2) = X( )(1,Q) which is of order 
r~ 2 , in the dipole case we start from £>(_3) = X( 1 )(/c 1 P 1 (z), Q) = X^-, f Q) + 
X(i)(u, §Q). The last two components are, respectively, of orders r~ 3 and 
r~ x . To compute the next terms we use operations curl" 1 and While 
curl" 1 generically increases the order by one, ( l39i) -( )40l) show that # applied 
to X( n )(/, Q) or Y( n )(/, Q) leads to terms of different orders: in generic case, 
first three terms of (1391) or (HUj) have order increased by one, the next two 
terms there have order unchanged and the last term in (HOl) has order in- 
creased by two. We present the results of computation in the monopole and 
in the dipole cases in terms of u = ar/2 and R = log(j^) = —2u + . . ., 
obtained using the symbolic calculus provided by the program MAPLE 8. 
The terms of order higher than r (or u) are denoted by o(u). We collect all 
the terms of generic order r l from different parts of fH3|) or f|44|) and denote 
them, respectively, by T>^ or B®. In the monopole case we assume that the 
charge Q = const. 



Q 




ii 



■2n+l 



The monopole case 



p = X>(- 2 ) + +o(u), B = £ (0) + o(u) 



(49) 
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where 



P(- 2 ) = X (0) (1,Q), D (1> =X (1) (/„iQb), B^ =Y (1) (u,-Qb), (50) 

cr a 

with 



f = --{u + u- v )R 2 -R-u = -u 3 + 



(51) 



The dipole case 



where 



V = _^Q 5 (3) + p(-3) + P (-i) + P (o) + P (i) + o(m)) (52) 
3 

B = B ( -V + + + £J (1) + o(«), (53) 



2>(- 3 >=X (1) (i*-\|q), (54) 
V™ = -X (1) (A, + X (1) (u, |Q), (55) 



^ (0) = X (2) (/ 2 , ^(b V Q) ) + Y (1) («, ^(b xQ)') + 



u 



+ Y (1) ( — ,lbxQ) + X (2)(/3 ,^bVQ), 



(56) 



= - Y (2)(^' ^ b >< ( b V Q)) " X 0)(/4, V (b V Q)) 



-Y( 2) ( 



u 



2a 2 



bV(bxQ))-X (1) (/ 5 ,-Q) + 



-X (1) (/ , ^(aQ + |Q)) - X (1) (/ , ^(dQ + |Q)' ) + X (1) (/ 8 , ^Q) + 



X 



(i) 



d Q 1 dQ 



+ X (1) (/ 6 ,J,(g ) " ) ,(57 ) 



B (-D 



b<- 2 > = -y (1) («-\|q), 

^Y (2) (l,bVQ)+X (1) (u,^bxQ), 



(58) 
(59) 
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B<°> = -Y {1) (u, aQ + |Q) + Y {1) (/ 1; A" ) + Y (1) (/ 9 , -^Q), 



= ^Y (2) ( M 2 , b V Q) - Y (2) (/ 10 , ^(bVQ)')- 
Yp,(/ 2 , [^(b V Q)" J" ) - X (1) (/ , ^j(b x Q)- ) - 
X(i)(/o, i[^(b x Q) ]• ) - Y (2) (/ n , £jb V Q) - 

Y (a) (/s, (^b V Q) ) - X (1) (/ 12 , ^bxQ)- 

X(1)(/5 ' ^ ( ^ b X } " Y(2)(/l3 ' ^ b V ^ + 



X (l)(/l4>^ b X Q), 



where / is given by fl5Tj) . 



1 



/j = --(m" 1 + m)# 2 - R + u = 2m + . . . 



h 



16 V ' 4 



i? 2 --(n- 1 + M )i?-- = M 2 + ... 

4 V 7 4 



/ 3 = --(3m- 2 + 2 + 3M 2 )i? 2 - -(u- 1 + - 3 5m 2 + 6m 4 = 2 
^ 8 V ; 2 V ; 2(1-m 2 ) 

/ 4 = - — {5u~ 3 + 3u~ 1 + 3u + 5u 3 )R 2 --{15u- 2 + U + 15u 2 )R 
32 8 

15U" 1 + Au - 23m 3 + 20m 5 _ 10 3 

8(1 — m 2 ) + '"' 



/ 5 



M 



1 — W 



-U 6 + . . . , 



/ 6 = ^(m" 1 + u)R A + -R 3 - -{u- 1 + M)i? 2 — R — u — -2m 3 + 
48 6 4 



./V 



1 _a 2 2m + 4m 3 

-Ti -_ri — 

2m 1 — m 2 1 — m 2 



-2u 3 + 



fs = ^{u- 1 + u)i? 4 + ii? 3 - ^{u- 1 + u)R 2 - 
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3 - 4m 2 + 3m 4 „ 3m - 3m 3 + 2u 



5 



R ^ = ~2u 6 + ... 



2(1 -W 2 ) 2 2(1 ~« 2 ^ 

.-1 .Ad2 



1 . , , 9 1 + u 2 3u — u 3 
4 v 1 - u 2 1-u 2 



3 



Ji3 = -7(1 + 3u z )R z - —R H - ? - = 6m z + . . . , 



4 1-u 2 1-u 2 



7 Description of the functions appearing in 
the Maxwell fields (a conjecture) 

The method of analyzing Maxwell equations in a neighbourhood of a freely 
moving electric multipole leads, in case of a monopole and a dipole par- 
ticle, to an interesting family of relatively simple functions of the variable 
r. These functions are obtained from A n (r) = P n (z) contained in (I46j) 
by consecutive application of the differential operators contained in # and 
the integral operator contained in curl" 1 . The latter is based on solving 
the differential equation fl30l) and could a priori lead very quickly to non- 
elementary functions. An unexpected result of our analysis consists in the 
fact, that the functions arising here are polynomials of the universal quantity 
R = log(j^) = tog(|qp^), with coefficients being rational combinations of 
u = ar/2 itself. This fact leads us to a conjecture, that a new family of 
special functions arises here, which makes the physical picture of the field 
relatively simple. 

In this Section we present a conjecture about the smallest algebra of func- 
tions / appearing in X( fc )(/, Q), Y^(f, Q), when we solve Maxwell equations 
with multipole particles, using the method proposed in this paper. 
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We set R s = First we notice that only constant functions appear in 
X(o) (remember that Ym) = 0). Indeed, the leading term for the monopole 
particle contains only Pq{z) = 1 and the terms with X( ) produced by fl39|) - 
( 110|) vanish. 

We define H k as the linear span of functions / entering ~X.( k )(f, Q), Y^(f, Q), 
which appear when we consecutively use formulae ( |32l) - ( |33j) or ( 1391) - ( HOI , 
starting from multipole fields ( )46l) . More precisely, we take the following 

Definition. We set H k , k = 1, 2, . . ., as the smallest vector spaces such 
that 

u- k ,u- k+2 ,...u k e H k , (62) 
and such that for any / G H k one has S(fc)(/) £ -f^fc, 



1 -m 2 



g if fc , g fl- fc> (l - u z ) [u z f, uu + 2uf >u - k(k + i)f] g H k , (63) 



2u 3 f 
l-u 

2u 3 f 
l-u 

Remark. Due to condition (1621) we have: 



u 2 f, u + (k + l)uf G H k+1 , u 2 f, u + (k + 3)uf + J - G H k+1 , (64) 

- u 2 U + kuf G if fc _ ls u 2 /,« - ik - 2)uf + e Hk ^ (k ^ 1). (65) 



P n (z) = Pni^iu- 1 + u)) G span{ M -", u~ n+2 , ...,u n }, 

and, whence, functions / contained in £>(_ n _ 2 ) = ~K^(k n P n (z), Q) are in H n . 
Also for n = and / appearing in = Y(i)(w, ^Qb) we have f = u £ H\. 
It remains to consider operations curl" 1 and # which yield B(i) and P^) with 
higher /. Due to (j52} - (j33jl and ([MD (SOD, conditions (E3D (ESD imply that 
also / appearing in those Bm, Vn\ belong to the corresponding H k . 

It is easy to see that conditions (IB"3"1) - (I5"51) can be reformulated in a simpler, 
equivalent form: 

1 G H k , uf )U G H k , (1 - u 2 )[u 2 f, uu + 2uf )U - k{k + l)f] G H k , (66) 
ii f 

" eH k+u u 2 U + (k + l)uf eH k+ll (67) 



1 — u 



2 



l-u 
uf 

eH k ^u 2 f jU -kufeH k ^ (k^l). (68) 



l-u 2 
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Conjecture: The space H k can be represented as 

H k = I^J H km (increasing sequence of vector spaces), (69) 



where H km are the direct sums: 



Hkm = H k j m , (70) 

3=0 



^ = ^((^^'(^^'■■■'(^31)^;' m = 0,1,2,..., 

(71) 

f R 2l-l u k+l R 2l-l u k+3 R 2l-l u 2m-k-X} 

H k , 2l ^ m = span (T^rp TyT^Tp - , (|f2 _ 1)m ) , (m > k), 



(72) 



-^fc,2«,m = span 



w l k (R,u) w l k _ 1 (R,u) w[(R,u) 



u k [u 2 - l) m ' u k - 2 {u 2 - l) m ' " ' ' u- k+2 {u 2 - l) m 

R 2l u k R 2l u k+2 tflrfrn-k w [( R , u -l) u 2 m -k+2 



(u 2 - l) m ' (u 2 - l) m ' " ' ' (« 2 - l) m ' (w 2 - l) m 

w l 2 [R,u- 1 )u 2m - k+4 w l k (R,u- l )u 2m+k 



(u 2 -l) m ' ' (w 2 - l) m 
Z = 1,2,..., with w l k given recursively by w l (R, u) = R21, 



(m > k), (73) 



(rain(fc,2Z-l) \ 
J2 A kj w[zf 2] (R,u)u 2 i- 2 \ +(2k-3)\\R 2l ^u 2k -\ (74) 

k, I = 1,2,... (here, we put (—1)!! = 1, so that (—1)!! • 1 = 1!!), for some 
real numbers A k j, k = 1, 2, . . ., j = 1, 2, . . . , k. Moreover, we conjecture that 
A k \ = 2k — 1, A k2 = 1, A kk = (2k — 3)!!. Unfortunately, we were not able to 
find a general formula for A k j. 

Remark 1. We know only a recursive procedure to find A k f Suppose 
we have found all A r j, r < k. Then we can compute w l r , r < k (cf. (J74"|) ). 
We set 

_ w l k _ 1 (R,u- 1 )u 2m+k - 2 

f ~ ( M 2 _ !)m ' 
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i.e. / is the last but one element of (175|) . Then we compute g G H k by means 
of the last formula of (JHSJ). Next we decompose g into the bases of H k j tTn+1 , 
j < 21, where w l k is given by f!74|) . This gives us an equation which may 
be used to find coefficients A k j, j = 1,2, ... ,k. That kind of procedure was 
implemented by us for k < 10, using the symbolic calculus provided by the 
program MAPLE 8. In particular, we have obtained the following results for 
k < 6: 

A n = 1, (75) 

A 21 = 3, A 22 = 1, (76) 

A 31 = 5, A 32 = 1, A 33 = 3, (77) 

A 41 = 7, A 42 = 1, A 43 = 50, A 44 = 15, (78) 

A 51 = 9, A 52 = 1, A 53 = 1273/9, A 54 = 485/9, A 55 = 105, (79) 

A 6 i = 11, A 62 = 1, A 63 = 62564/225, A m = 2703/25, A 65 = 9985/3, A m = 945. 

(80) 

Using (174|) . one gets that 

min(fc,2«-l) 

wi(i2,«)= 9ks{u)R 2 i- s , (81) 

where ^fc s are polynomials. In particular, we have: 

w l (R,u) = R 2l , (82) 

w[(R,u) = R 2l + uR 2l _ l , (83) 

U4(#, U) = 3R 2 l + (3M + M 3 )i?2i-1 + ^ 2 i?2i-2, (84) 

w l 3 (R,u) = 15R 2l + (l5u + 5u 3 + 3u 5 )R 2l _ l + (6u 2 + 3u i )R 2l „ 2 +u 3 R 2l _ 3 . (85) 

(Observe that, due to (IHTT) . for / = 1 we omit the last term in w l 2 and the last 
two terms in w l 3 .) 

Remark 2. One can prove that for A^i = 2k — 1 and any value of 
the remaining coefficients Ay, the elements (17Tl) - fl73l) . spanning each vector 
space Hkjm are linearly independent and the direct sum condition is obvious. 
Moreover, using (1741) . Hk m form an increasing sequence of vector spaces and 
therefore Hk are well defined. The elements f )62l) belong to C Hko m C 
Hk (m > k). Using once again f!74|) . it can be also checked that the first 
operations of fl6"6"l) - fl6"8l provide elements of H k , H k+ i and H k _i, respectively. 
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It remains to find the formula for A k j such that the remaining four relations 
of (Kffi|) - (K)8|) are fulfilled, to prove that s^)(f) G H k for / e H k and to show 
the minimality of H k . 

Remark 3. Elements of H k are rational functions of u and R. However, 
the generators of these spaces, proposed in formulae f!7T]) - fl73|) are quite com- 
plicated. But splitting the nominators into monomials doesn't simplify the 
situation, because application of s^) (which must preserve H k ) to the com- 
ponents obtained in that way leads to functions which seem to be in general 
not rational in u and R. Thus the form (!7T|) - (!73|) of the generators seems to 
be the simplest one. 



Appendices 

A Properties of fields X and Y. 

We assume r 7^ 0. Using (|2"5I) - (12"11) . after some computations we get 
d t X k = i-r 2 l„ + (2n + 3)rf >r - (n + l)(n + 3)/]^§tf+ 

[(» + !)/ - rf,r)^6 u H + n[(n + l)f - rf >r ]^+ 
(rf, rr - nf, r )^§ + ^(n - l)Q kh 



r r fr a 2 r 2 f (n + l)f,xiE k 
diY k = + — -f - K - - 44+ 

-Fkh 



where H Qii...i n -£ x . . . x , Q k Q ki2...i n x ■ ■ ■ x , Q k i Qkii^...i n x 3 ■ ■ ■ x 
E k = e k ' m XjQ m , F k , = ( k ' m XjQ m ,. That implies (2SI) (Z 

Xi x>=<2±Mh, x«y' = o 

and (cf. (TT6J)) 

X k=2l r Ur + rf,r -{71+1} /] + 



r n + l - - (86) 
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-^[rf,r + {n + - 2 [(" - 2 ) r /,r ~ ™(™ + 1)/] rn+1 + 



-[rf, rr + (n + 2)/, r ] ^ + 2 rn _ 2 



l rU oV/ (w + 3)/ g fc (6*s t ) 
V </? 2<^> 2 J r™ +1 

,n-l K O/nr-n-1 KL mr-n+l V 1 



where G k = e k lm biQ m . Let us notice that Q m m = (Q is traceless) implies 

= 25 ms Q ms e al k x a b l = D ijs e ijs e al k x a b b 
where D^g = €ij m Q ms . That and 

tijstaik = 8ia(8ji5 S k — Sjk6 s i) + 5u(6jkS sa — Sj a 5 sk ) + Si k (5j a 6 s i — 8ji5 sa ) 

give 

e alm x a kQ mk = -b l F kl + G k . (89) 
Using and fl89]), we obtain 

X k(n+1) (a, bVQ) = [(n + 2)a - ra^K&'x^ - ^-^r^QJ] + 



(n + l)(2n + l) 

n (n )(/3,bxQ) = ^ I ^[n6 fc ^-(n-l)x fe (g m 6 m )-n(^)g fe +(n-l)r 2 6'g H ], 
rvpr n+L 

(91) 

^(n-DC^bJQ) = (w7-n, r ) Xfc( r 6 ^ m) +^hb m Q mk (n>2), (92) 

n ( n + i)(^,bVQ) = ^ 2 [-^-(e k %b l )H+-^-E k (b l x l )- - " ( " ~ \- r% l F kl \ 
tpr n+z n + 1 n + 1 [n + l)[2n + 1) 

(93) 

X k{n) (K, b x Q) = ^[G, - ^b l F kl ] - [(n + 1)k - r «, P ]^§£ (94) 
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Y k(n „ 1) (p,h\Q) = -^b l F kl (n>2). (95) 
Computing e km e m li ' ej ac x k XiQ a b c by means of ei k m e mij = S H S kj - 5ij5 ki or 

£mi J £jac = S ma 5 ic ~ $mAa, One obtains 

(e^Xjb^H = (b [ Xl )E k - x k (b l E{) - r 2 G k . (96) 

Proof of Theorem 1. Comparing (l87]) - (l88l with (!90|) - (l95|) and using (|96|) . 

we get 

X k = X k{n+1) (a, b V Q) + Y k(n 0, b x Q) + X k{n _ 1} ( 7 , bj Q) , (97) 

Y k = Y k[n+1) {ijj, b V Q) + X k(n) {K, b x Q) + n (n _i)(p, bJQ), (98) 
for a = -\[r 2 f, r + (n + l)rf], (3 = g 1} 7 = 2r £~l 1) (-r 2 f, r + nrf), if; = g 2 , 

1 fr 2 
K = -2^+1) V' P = 93 ' 

Using (|23p - ()24p . we obtain that X k gives the last two terms in fl34l) while 
Y k gives the last three terms in flM}. That, (TJ5]) and ([9"7]) - (j9"8"ll prove (153)1 - 
dMD. Q.E.D. 



B Proof of the properties of A n and 5 n 

The quantity 

1 ( 1 ' ' 2 k ~ l a- k r~ k 4 



2 \ ar 2a 



is of order r~ k . Therefore P n (z) (with the highest term of order z n ) is a 
Laurent polynomial of order r~ n and analogously t> n _i is of order r~ n+1 . 
Moreover, 

1 , z + 1 , 1 + w 



log ■ 

2 z 



log 



1 



ar + 



(99) 



and, therefore, B n (r) can be written as a Laurent series in r, starting from (at 
least) r~ n+1 . Inserting this series into h^(B) = 0, using (1281) and denoting 
the order of B n (r) by Z > —-a + 1, we get that vanishing of the r l ~ 2 term in 
h( n )(B) implies I = —n (which is impossible) or I — n + 1. Hence, B n (r) is 
of order r n+1 . 
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C The transformation of 2 n — poles (n = 0, 1, 2) 

Let X A be the current density in the Fermi frame (coordinates £) and J7" M be 
the corresponding current density in the modified Fermi frame (coordinates 
x). Then 



J 11 - 



det 



dx 



Integrating with the test function f(x), we obtain 

S" = J J»(x)f(x)d 3 x = j^.X^\det(^)\f(x)d 3 x = J ^Z\t)~f{t)d% 

(100) 

where /(£) = f(x(£)). Differentiating ((31), one obtains 

_ = (Sf + a%)M - (e + -a k p 2 )( ai + -a%)M\ (101) 

a k p 2 M - {£ k + -a k p 2 ){a m C + -a m a m p 2 )M 2 , (102) 



d£° 2 r ^ 2 r /v 2 

ia 2 p 2 ) -1 - Moreover, g£ = 0, w 



where M = (1 + a k £ k + Wp 2 )' 1 . Moreover, f£ = 0, g£ = 1. Thus 



5° = | T°(0f(0d% (103) 

where 

^ - f ^ + J*- d05) 

Suppose that X A = 2^\(£, Q) is the standard 2 n -pole, given by the right 
hand sides of the formulae (|4ip - (j42p with respect to the coordinates We 
shall calculate the corresponding value of = J~f\ by means of fllOOp in 
the case of a monopole, dipole and quadrupole, i. e. for n — 0, 1, 2 (for n = 
we assume Q = const.). 

For n = T ( ° 0) (£,Q) = AtcQS^, Jf 0) (f,Q) = 0, hence S° = 4ttQ/(0) = 
4ttQ/(0), L fc = 0, S fc = 0, 

Jg ) =25 ) (a:,g), (106) 
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where T^Jx, Q) is given by (jH|) - (Jl2"|) with respect to the coordinates x M . 

For n = 1 I° (1) (^Q) = -AnQ k d k 5^\ Jf^&Q) = ^Q k 5^\ hence 5° = 
A7iQ k (d k f)(0) = A7tQ k (d k f)(0), L k = 47iQ k 5^, S k = 47rQ fc /(0) = A7tQ k f(0), 

J^=l^(x,Q). (107) 

Forn = 2 2f 2) (£,Q) = f Q^d^S^, 2f 2) (£,Q) = -^fQ kl d^\ hence 
<?o = ^ Q M {dkd j m = ^ Q ki {dkdifm _ ^QW a ,(9*/)(0), L fc = 
-f Q fcm a m 5 (3) - f Q fcm a m <5( 3 ) - ^Q«Si<j(3), S fc = -M(Q« aj )- /(o) + 

£q*W)(o), 

^ ) =Jf 2) (x,Q)+Jf 1) (a;,P), (108) 

where the dipole charge P fc = — ^Q kl ai. Therefore the Fermi frame quadrupole 
has a dipole component in the modified Fermi frame. 



D Integrals, distributions and the proof of 
Theorem 2. 

First we compute the following integral over two-sphere: 

S kim = J ) (^) ^ (f) 21 (^) ^ da, k, I, m — 0,1,2, ... , (109) 

where r = (x 2 + y 2 + z 2 ) 1 / 2 is the radius of the two-sphere, da = sin 9d8d(p, 
r, 9, 4> are the spherical coordinates. 
Proposition 1. One has 

(2k- 1)!!(2Z- l)!!(2m- 1)!! 
klm ~ (2A; + 2/ + 2m + l)!! 

Proo/. Let A > 0. We use 

P POO 

/ x 2k e- Ax2 y 2l e- A y 2 z 2m e- Az2 dxdydz = / r^+^e-^V^dr, 

JR3 JO 



x 2k e~ Ax2 dx = (-d/dA) k [ e' Ax2 dx = ( 2fc ^1/2^-^-1/2 
'r Jr 2 

and the analogous formulae for y, z and r. Q.E.D. 
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Corollary. 

r Airr 2n . 

Proof. We get nonzero results only if i s = 1 in 2k cases, i s = 2 in 2Z cases, 
z s = 3 in 2m cases, + I + m = n. Then the left hand side gives Sum r2n 
and the sum gives the factor (2k — 1)!!(2Z — l)!!(2m — l)!!2 n n!. Next we use 
(HDSD. Q.E.D. 

Proposition 2. Let Q ll - %n be a symmetric traceless tensor, S^"^ m be 
any tensor. We set Qx ■ ■ ■ x = Q %1 "' ln Xi l ■ ■ ■ Xi n and similarly for S. Then for 

m = n 



(Qx ...x)(Sx---x)da = Q^ Sll ... tn -^—r 2n , (111) 
S 2 (r) (2n + l)\ 

while for m < n the left hand side of (111 II) equals zero. 

Proof. Let m = n. Then the left hand side of fillip is equal to Q l ^- l ^S ln+1 - %2n 
multiplied by (jllOp and summed over all We may first sum over all it 
and then over A € IT^n- In such due to the traceless condition for Q, 

nonzero terms in the sum over A are obtained if for each k = 1, . . . , n one of 
elements X(2k — 1), \(2k) belongs to {1, ... , n} and the other one belongs to 
{n + 1, . . . , 2n}. We get 

A— Jin 

O n n\ D* 1 -*™ Qh-3n X. X. 

(2n+lV Z^Z^ V %(1)J1 ■ • -%(n)Jn) 

which due to the symmetry of Q gives the right hand side of (II lip . For 
m < n a similar arguments show that all terms vanish (we use the traceless 
condition for Q if m+n is even and the antisymmetry of the expression under 
the integral if m + n is odd). Q.E.D. 

Let n = 0,1,2,. .., Q be a tracelless symmetric tensor of rank n, f be an 
even analytic function of r near (/ is a constant for n = 0). We set 

„ „ Qx ■ ■ ■ x 

n Q/ - / r .2n+l 

and define T as the linear span of all -DnQ/- The elements F £ J 7 become 
distributions if for any test function p we set 



oo 



<F,p>= dr r 2 Fpda = lim / Fpd 6 x. (112) 

J0 JS 2 (r) R ^ 0+ Jr 3 \K{0,R) 
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That is well defined because setting F = -D„q/ and using 



n— 1 -. 

Ph...i k = (d h ...d ik p)(0), 



we obtain 

r 2 Fpda 

S 2 {r) 



n-l 



/ j fa\ y2n—l 



k=0 



(Qx ■ ■ ■ x)(px ■ ■ ■ x)da 



S 2 {r) 



(113) 
(114) 

+ 0(r)=0(r), 



since all the integrals over S 2 (r) vanish due to k < n and Proposition 2. 

For F G T let <9j denote the partial derivative in the sense of distributions, 
df - the partial derivative as function, 



df-F = diF - d?F. 

We define 1; as tensor of rank 1 such that = 5ij. 
Proposition 3. 

Tl 

XiD n Qf = D n+lt i iV Q tr 2f + —D n _ lt i^Qj, 



;ii5) 



2n + 1 



@i D n Qf — -D n +l,l I VQ,-(2n+l)/+r/,', 



n 



2n il^-iMQilr/r- 

Proof. We use Xi(Qx---x) = (lj V Q)(x ■ • • x) + -^^Qi-, where Qi = 
Q ii2 .., n x i2 ...x i " = (l i \Q)(x...x). " ! Q.E.D. 

Thus for F G T , df F and df'F are distributions. 

Proposition 4. For F G T and a test function p 



< d*F,p>= lim 



FpXiR da. 



Proof. Due to (fTT2l 

< diF, p >= - < F, dip >= lim 



(dfF)pd i x- lim 



R ^ 0+ 7R3\ir(0,R) 



df{Fp)d i x, 



R ^ 0+ JR3\K(0,R) 



where the first term equals < df ' F, p >. Setting {Gi)j = Fp5ij and using 
(11151) . one gets 



< d*F,p>= - lim 



divG,- d x 



R ^ 0+ JH s \K(0,R) 
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= lim I (G i ) i —R 2 da= lim / FpXiRda. Q.E.D. 
R^+J S 2 {R) R R-*o+Js2 {R) 

Proposition 5. 

d«D nQf = fe/(0](-l) n - 1 ^y |j ^5( 3 ), 

where I = (i 2 , . . .i n ), dj = d ia . . . d in . 

Proof. Using Proposition 4 for F = D n Qj, (11131) . setting c = lim r ^ fi r ), 
^p.yi...i k i k+1 _ pti-.-ikfii-k+i an( j us j n g Proposition 2, we get 

Pii% ' ' ' P %ii ' ' ' •Kik^i] 



c _ iia 4 47r(n!) 2 2" 
;Q ,ta '" tn Pi a ...i n - 



(n-l)P rt2 - In (2n + l)! 

Moreover, we use < d/tf^p >= (-l) n_1 Pi 2 ...i„ (cf. (HI) and ^g+ ! = 
(2n + l)!! " Q.E.D, 

Proposition 6. Assuming / G L n0 (cf. Section 4) and setting g = fr n 
(which is regular at 0), we obtain 

71+1 

-Xfc(n)(/)Q) — Ai+l,lfcVQ,(2n+l)s-T-g, r + 2 ra _|_ \ •^ n - 1 > 1 feJQ,fl,r/n 
Yk(n)(f,Q) = -Dn,l k xQ,g/<fi- 

Proof. By a direct computation. Q.E.D. 



Proposition 7. 

9fX fc(n) (r-,Q) = (-l) 



r v , -„ ^ , lV1 47r(2n + l) 



(2n + 3)!! 



d«X k(n) (r-^\Q) = (-^T^^Qik L d L S^, 
d?X Kn) (f, Q) = if hm(/r"- 2 ) = 0, 
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dfn (n) (r-",Q) = (-If — - - [(n - l)e k ma Q mi L d aL 5® + e^Qjdjd^}, 

Q)=0 if lim(/r») = 0. 

r— »0 

Proof. We use Proposition 6 and Proposition 5. Q.E.D. 

Remark. In the case of n = 0,1 some equations related to the first 
equation of Proposition 7 were presented in (6)-(7) of 

Next, we obtain that in addition to ([25]), Q2EH2B, (|HH35P, the op- 
erations div, curl^ and # acting on the fields X and Y yield the following 
distribution parts: 

Proposition 8. 

div«X M (r-»,Q) = (-l)"|^±ilQ^P), 
[curlJX (n) (r-,Q)] fe = (-1)" + e^Q^M, 

[curl*Y w (r-,Q)] ft = (-1)"^^^^ 

div B Y (n) (r- n , Q) = Y (n) (r-", Q) #/? = 0, 

while for / such that lim r ^ (/ r?1 ) = the distribution parts vanish. 

Proof. We use Proposition 7 and notice that the operator d/dr appearing 
in # (cf. (|15[) ) gives no additional terms. Q.E.D. 

Let us notice that tpo = -Don = - is the classical monopole potential, 
-A^o = 4tt(5 (3) . By induction one has Q 7 (9/^o = (-l) n (2n - l)H£> nQ1 (no 
distribution parts! - cf. Proposition 5) and therefore i[) n = D u qi satisfies 

4-7T 

" A * = (- 1 )"(2^ <^ '* ^,3, 

(cf. OH])) and ip n is the classical 2 n -pole potential. Using ( 11 15ft and Propo- 
sitions 3,5 and 6, we obtain that the corresponding electric field is given 
by 

l£ = -grad^n, (V%) k = X k(n) (r- n , Q) + (V d ) k , (116) 
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where 

^»' = (- 1 »"(2^ W3) 
(cf. (145p and Remark 2 in Section 5). After some computations one obtains 

^^(-ir^^Q^K (H7) 

(curl^)* = (-l)^ 1 ^ 4 ^)!! ^ 3 '^"^' ( 118 ) 

= (" 1 ) n ?2^ ! T)!!^ 5(3) - ^ 

Proof of Theorem 2. Let us first omit the terms containing 5^ and its 
derivatives. Then Dm and Br a consist of X and Y fields and therefore (cf. 
fl25|) ) divTfy) = divi3(j) = 0, which proves f|T7|) - f|T8|) . Moreover, applying 
curl^ to (J4TD (J48D, we obtain 

curl v B( k+ x) = k=-n-2, -n, -n + 2,..., (120) 

- curyD (fc+1 ) = i3 (fc ) # , k = -n - 1, —n + 1, —n + 3, — (121) 
But one has also 

curyD ( _ n _ 2) = (122) 

(we use (T27|) and ft, (n) [P„(z)] = 0). Inserting fjl3"]) - (13I]) into the both sides of 
(JUD d2QD and using (TT20]) -ffT22l). one proves that (pi)- (f20l are fulfilled. 

Let us notice that 5^ and its derivatives appear in the sources (14T1) - (142|) . 
in D d and also when acting by div, curl v or # on X( fc )(r _fc , P), Y( fc )(r~ fc , P). 
Using (I32l - (l35l) . one shows that X( fc )(r _fe ,P), Y( fe )(r~ fe ,P) appear only in 
P(__ n _2) which contains X( n )(r~ n ,Q) and in S(_ n _i) which contains 

- £tt Y(n - i)(r_(r " 1) ' bJ Q) " Y w (r_n ' ^ 

Therefore div£> contains additionally 

4-7T 

divD, + div fi X n (r-™, Q) = (-1)"- 7^Q J dj6® = J°, 

(Zn — 1)11 
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div£> has no additional terms, X># — cml^E contains additionally (cf. Propo- 
sition 8, (jUSJ) and (S2D) 



P/ + X (n) (r-",Q) # « + 



n 2 - 1 
2n + 1 



curl^Y^Cr-^.bJQ) 



+curl£Y (n) (r 



— n 



Q) 



B* + curl^P contains additionally curl^D^ + curl^X( n )(r~ n , Q) = 0, hence 
the Maxwell equations are satisfied also in the distributional sense. The last 
statement of Theorem 2 follows from the properties of curl" 1 and X and 



We thank Dr Marcin Koscielecki for fruitful discussions. 
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